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The daily global-mean values of the lower-tropospheric temperature determined from microwave emissions 
measured by satellites are examined in terms of their signal, noise, and signal-to-noise ratio. Daily and 30-dav 
av erage noise estimates are reduced by almost 50*3? and 35T. respectively, by analyzing and adjusting ( if necessary ) 
for errors due to 1 ) missing data, 2) residual harmonics of the annual cycle unique to particular satellites. 3) 
lack of filtering, and 4 ) spurious trends. After adjustments, the decadal trend of the lower-tropospheric global 
temperature from Januarv 1979 through February 1994 becomes -0.058°C. or about 0.03°C per decade cooler 
than previously calculated. 


1. Introduction 

Microwave Sounding Units (MSUs) onboard 
NO A A polar orbiting satellites have provided temper- 
atures of deep layers of the global atmosphere since 
1979 [Spencer and Christy ( 1992a,b, 1993 )( hereafter 
SC92a. SC92b. and SC93)]. These studies have dem- 
onstrated the high quality of the MSU data and their 
value in addressing global climate issues [see also Hur- 
rell and Trenberth ( 1992). Basist et al. (1995)]. In 
general, such studies have used MSU data averaged 
into monthly timescales. In this paper, we shall ex- 
amine the signal and noise characteristics of the MSU 
channel 2 retrieval of the lower-tropospheric temper- 
ature. know n as MSU 2R, for daily global- mean values. 
Our goal is to provide a dataset of daily global tem- 
peratures for which errors, or noise, have been reduced 
as much as possible while preserving the magnitude ot 
the signal variability. 

The time senes of daily zonal-mean temperatures 
through February' 1994 has utilized measurements 
from instruments on seven different satellites ( TIROS- 
N. XOAA-6. -7. -9, -10. -II. and -12). How these data 
were generated and merged into a single time series 
will be explained in the following two sections. To as- 
sess changes made to the original data, we shall begin 
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our analysis from the point where no quality control 
has been applied. We shall examine the impacts ol 
missing data, unique satellite characteristics, filtering, 
and spurious trends on the global-mean values and 
determine how to enhance the quality of the time senes. 
Users who have performed research on current editions 
of the MSU 2R data should realize that some quality 
control was applied to those data and that the final 
results of this study do not alter the earlier monthly 
versions in a significant manner. 

2. MSU temperatures 

Temperatures are produced from the raw MSU data 
in the following way. The radiance counts observed by 
the MSUs are converted into brightness temperatures 
by an interpolation scheme, which is anchored on the 
warm side by an onboard hot target ( measured by a 
set of redundant platinum resistance thermometers) 
and on the cold side by the temperature of cold space 
(2.7 K). On each scan, the sensor obtains radiance 
counts for these two targets of known temperature as 
well as eleven footprints over the earth. The counts 
observed in the earth views are converted to temper- 
atures using the interpolation scheme derived from the 
cold and hot targets (Spencer et al. 1990). 

As described in SC92b, a limb-retrieval technique is 
applied to the scan-line data, which results in a vertical 
temperature profile ot the lower troposphere. For this 
study, the daily data are binned into 2.5° latitude 
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Fig. 1. (Top) Data availabilitv Irom the operational NOAA polar orbiters, which earned Microwave Sounding Units. 
( Bottom) As above but categorized by crossing time: morning or afternoon. 


bands, which are then averaged into zonal-mean tem- 
peratures. Global-mean values are calculated from 
these zonal quantities. Comparisons with individual 
radiosonde locations indicate the MSU 2R monthly 
temperature anomalies correlate most highly with the 
1000-500- or 1000-400-hPa thickness anomalies 
(SC92b). The use of daily zonal means provides a ro- 
bust dataset from which to calculate global means rel- 
ative to daily gridded data, which have coarser reso- 
lution and much greater noise. 

The time period examined here is 1 January’ 1979- 
28 February’ 1994. For this 15+ year period, over 97% 
of the possible days had at least 85% of the possible 
observations recorded by at least one satellite. Most ot 
the days below the 85% daily data rate occurred be- 
tween 1979-81 (inclusive). 

3. Merging method 

A brief description of the method by which the data 
from seven satellites were merged is contained in the 
Appendix of SC92a. In the following discussion, we 
shall provide background information and an update 
of that method. 

There are basically two sun-svnehronous satellite 
orbits, w hich meet the operational requirements of the 
National Weather Service. One is the “terminator" or- 
bit in which the satellite monitors the earth at local 
time ( LT ) near sunnse and sunset, crossing the equator 
at about 0730 LT (descending or southbound) and 
1930 LT (ascending or northbound). A platform in 
this orbit is commonly called the “morning'’ satellite. 
The second or “afternoon” orbit crosses the equator 
at about 1 330 LT ( ascending ) and 0 1 30 LT ( descend- 
ing). The crossing times for trie morning satellites have 
been very stable through the years, while those ot the 


afternoon orbiters have tended to degrade, ascending 
later in the day as length of service increases. In par- 
ticular, the most recent platform, A r O,A.4-l 1, has drifted 
by 3 hours, creating potentially serious consequences 
as described later. 

Since 1979, there have been three reliable MSUs on 
the morning satellites: A ! O.A4-6, -10, and -12 The 
NOAA-8 data proved to be too noisy for necessary 
measurements, so the older morning satellite already 
in orbit (NOAA-6) was returned to operational activ- 
ities until A T OA4-10 was launched. None of the A O. L 1- 
8 measurements were used in any of our previous stud- 
ies nor will thev be used here. The lour atternoon sat- 
ellites are TIROS-N, NOAA-7. -9, and -11. 

The National Weather Service has generally 
launched a replacement satellite every 2 years or so. 
alternating between the morning and afternoon orbits. 
As such, the operational instruments onboard arc 
nominally anticipated to provide information over a 
4-vear period. This schedule has been interrupted due 
to such circumstances as A'O.TTS’s instability, a mal- 
functioning oscillator in the MSU on NOAA-9, and 
loss of an entire satellite, as in the August 1993 at- 
tempted replacement of NOAA-1 1 with NOAA-1 3, As 
a result, some instruments have been used for at least 
6 years ( NOA4-6 and - //). Figure 1 indicates the se- 
quence of launches and the data used in construction 
of the 15+ year MSU time senes (NOAA-8 included 
for display purposes). Note that there are some penods 
prior to 1988 in which only one MSU was providing 
data for the time series. 

The instruments are also arranged by crossing time 
category in Fig. 1. revealing the opportunities where 
simultaneous comparisons between MSUs on different 
satellites can be performed. It is in these overlapping 
periods where a robust measurement ot precision can 
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be made because the instruments on different satellites 
are completely independent. 

To combine the entire record into a single time series 
of daily anomalies, two forms of biases were calculated 
and then removed from each MSUs data in sequence. 
The first bias is due to differences in the mean annual 
cvcle as observed by the morning orbiters versus that 
viewed by the afternoon orbiters because they each 
sample the local diurnal cycle at different times. The 
greatest morning/afternoon satellite differences occur 
where surface emissions have large diurnal variations, 
since the total emission observed by the MSUs includes 
some emission from the surface. This is especially true 
over elevated terrain where the fraction of si:: ace 
emission becomes larger, as viewed from the satellite, 
because in such locations there is less atmospheric mass 
to produce emissions. Factors, which determine the 
amount of surface emission, are surface skin temper- 
ature and to a lesser extent, wetness, snow' cover, and 
vegetation. Variations in these parameters, particularly 
skin temperature, will disproportionately influence the 
MSU-observed brightness temperature over the course 
of a 24-hour day. 

In general, the ascending node of the afternoon sat- 
ellite (about 1330 local time) views sun-exposed land 
and the overlying atmosphere as being warmer than 
the other passes because of the increase in surface 
emission (warmer skin temperature) relative to at- 
mospheric emission. This is most apparent over dry 
regions with considerable sun exposure ( Sahara Desert ) 
or over high terrain (Tibetan Plateau) because day-to- 
night differences in surface skin temperature are much 
larger than those of the troposphere as a whole. This 
is a systematic effect and must be treated separately for 
the afternoon versus morning satellites. [The differ- 
ences in emission from water surfaces ( i.e.. oceans) 
between day and night is not a factor.] Therefore, to 
calculate anomalies of atmospheric temperature, we 
must determine the systematic biases of the afternoon 
MSUs versus those of the morning MSUs due to diur- 
nal effects such as surface emissions. 

We calculate a reference annual cycle based on Sep- 
tember 198 1 -August 1982, a period in which the 
morning ( A OAA-6 ) and afternoon ( NOAA-7) contain 
adequate data on nearly every day. Daily temperatures 
for this period are smoothed with a cosine filter (half- 
power frequency of 30 days) for each satellite, and the 
result is defined as the “preliminary" mean annual cy- 
cle. Because interannual variability exists, we are re- 
quired to utilize a period common to both an afternoon 
and morning satellite, which is at least 1 2 months long 
(i.e.. we can not use a different 12-month penod for 
the morning than for the afternoon satellite). Due to 
serious problems (to be discussed later). XOAA-I I is 
unavailable to serve as the afternoon satellite and only 
XOAA-7 contains a period of overlap longer than 1 
year with a morning satellite ( Fig. 1 ). Hence, the penod 
defined above is chosen for the preliminary mean an- 


nual cycle. As noted below, a longer base penod will 
be selected for the final product so that peculianties of 
September 198 1 -August 1982 will not unduly affect 
the data. 

We standardize all of the satellite data by subtracting 
the morning mean annual cycle from A OAA-6. - 10 , 
and -12 and the afternoon mean annual cycle from 
TIROS-N, NOAA-7, -9, and -//. In this way. we elim- 
inate the diurnal effects peculiar to morning or after- 
noon orbits from all satellite data. (The reader should 
keep in mind, however, that what has been removed 
are annual cycles as observed by only A OAA-6 and 
-7. These may differ in slight ways from the other in- 
struments had they operated dunng the 1981-82 pe- 
riod.) At this point, the total time series is a set of 
anomalies from the mean annual cycle as defined 
above. 

The second bias to be removed arises from the sys- 
tematic offset of one instrument versus another. 
Though the daily global anomalies of two simultaneous 
MSUs may correlate at 0.99+, each instrument has a 
relative temperature bias of as much as 1 K < S('92a. 
SC92b. SC93). We determine these systematic biases 
as the average difference between temperatures of the 
MSUs at each latitude during common periods of ob- 
servation. The bias is sequentially accumulated and 
removed for each MSU so that all anomalies then are 
relative to A OAA-6/7 for September I l - Vuc^st 
1982, our preliminary mean annual cycle. 

The final step is to select a representative jvnod as 
the base annual cycle and calculate anomalies reia! 
to this. We selected 1982-91 as the base because* * -e 
were several missing-day periods prior to and a 
length of 10 years seemed appropriate, though ot , ov 
it is arbitrary. We then determined the residuui an- jai 
cycle in those 10 years and subtracted this Mom no 
entire record. The result is 1 5+ years of daiK anomalies 
relative to the 1982-91 period. 

4. Precision during overlap periods 

Because the atmosphere is being observed h\ two 
independent instruments during the overlap periods, 
we may determine the signal and noise characteristics 
of the instruments. In Fig. 2 we show results of the 
mean and mean difference of the global anomalies ot 
the XOAA-I 1 and -12 overlap (hereafter 11/12). We 
shall define noise (A r ) as the variance of the anomaly 
difference divided by 2, N = (cr A/ 2 )“, an< ^ signal 
(5) as the variance of the anomaly sum divided by 2. 
S - ( o'-/:) 2 . The daily X of 11/12 is 0.005565 K~, 
while the" daily S is 0.032184 K 2 , giving a signal-to- 
noise ratio (57 AT of only 5.8. The noise, signal, and 
S/X for 30-day averages are 0.001211 0.016512 

K 2 , and 1 3.6, respectively. These are the lowest (worst ) 
S/X values of the six overlaps that are used in merging 
the dataset into a single time series. Improvements in 
5/A r will be most evident in this example. 
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Fig. 2. (Top) The average of the daily global-mean anomalies for 
SOAA-ll and -12 (Bottom) The difference of the global-mean daily 
anomalies (offset by — 0.8°C for clarity) divided by 2. This was the 
worst case ot the six overlap comparisons. 


Our goal is to reduce the noise without compromis- 
ing or smoothing the signal. We note that the MSU 
datasets previously provided to the scientific commu- 
nity contain some quality control aspects and, in real- 
ity, the final result of this present study will not ma- 
terially alter any conclusions one may have deduced 
using older versions of the MSU 2R data. We are 
showing, here, the results of data analysis for which no 
quality control was performed, and then we shall in- 
troduce measures that will optimize the S / N for both 
daily and 30-day mean statistics. The effects, which 
relate to noise reduction are 1 ) variations in the amount 
of missing data allowed, 2) mean annual cycle differ- 
ences between A OAA-6 or -7 and the other satellites, 
3) median latitudinal filtering, and 4) removal of 
NOAA-U's spurious warming trend. 

a. Impact of data availability 

The information for the statistics presented above 
for 11/12 were derived from MSU 2R data for which 
at least 20 000 observations per day were available from 
both satellites. For this section, we use as our “number 
of observations per day,’’ the total number of footprint 
observations from which the MSU 2R retrieval is de- 
rived. There are about 3000 useful scan lines per day, 
each of which provides eleven footprint observations. 
The MSU 2R retrieval scheme utilizes only eight of 
these observations for the scan-line retrieval, hence, 
the maximum number of distinct observations on a 
given day is about 24 000 (SC92b). The 20 000 value 
then represents about 85% of the total number of ob- 
servations that are possible. 

Since one complete orbit provides 7% of a day s data, 
errors are introduced by subsampling the earth using 
daily data amounts that represent fewer than 14 orbits. 
This is especially true where the surface emissions are 
strong. These areas require systematic sampling each 
day so as not to introduce biases into the global mean. 
Even missing only that part ot the orbit which views 
these “hot spots” can lead to greatly increased noise. 


We first attempt to reduce N by assessing how- miss- 
ing data impact global-mean temperature values. We 
have calculated VA' = ( 0 ^/ 2 ) l° r 19 000. 20 000, 
2 1 000, 22 000, and 23 000 per day data rates for daily 
variations and 30-day average variations ( Fig. 3 ).JThe 
values plotted are those of a weighted mean of VA' for 
all six overlaps based on the number of days in each 
overlap period ( recall we have mentioned the statistics 
of the worst overlap example earlier). The total number 
of days on which both of the intercompared MSUs 
reached acceptable data rates were 3009, 2895. 2697, 
2 1 33, and 1436 for 19 000 through 23 000. respectively. 

The results of the daily comparisons indicate that 
the noise is reduced as the number of observations re- 
quired for comparison is increased, as would be ex- 
pected. The reduction is noticeable but relatively slight 
with daily V/V decreasing from 0.068 K to 0.064 k tor 
19 000-22 TJOO data rates. There seems to he no sig- 
nificant difference in l/A f due to data availahilitv lor 
30-day averages, all being about 0.026 k Of all the 
characteristics we shall examine in this Mudv data 
availability has the weakest impact on A i '• m -us 
result, there is no compelling reason to hm • c k- 
ceptable data rate beyond 20 000, which u* 

currently used in monthly reports (e.g.. 
mate Diagnostics Bulletin). 

b. Residual annual cycle differences 

We have indicated that the morning a no u 
instruments observe the earth under diffenm - 
conditions, leading to different mean annua 
was assumed that all morning satellites wm..: , 

identical mean annual variations observed 
6. and all afternoon satellites those of NO. \ 1 
examination of Fig. 2 reveals a series ot v'e. . . 
in the difference time series, which happen t ■ • * • 
locked with the annual period. Because Fig - : * - 

the last overlap period, in which the MSI : r * * 
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Fig. 3. Solid (open) columns represent W for daily < 30-dav mean) 
statistics. The left-to-right progression indicates results tor minimum 
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Fig. 5. Solid (open) columns represent V ATor daily (30-da> mean) 
Fig. 4. Solid (open) columns represent VATor daily (30-dav mean) statistics due to median filtering with window width In + 1 indicated 
statistics. The left to right progression indicates results for removal t» e | ow 
of Fourier components of the annual period up to the wavenumber 
indicated below the columns. 

c. Median filtering 

adjusted with the accumulations of all previous ad- yh e nort h ern rnidlatitudes tend to have a large in- 
justments, any systematic errors should be most visible f| uence on ^e global-mean anomaly because of the 
here. Inspection of all overlap examples indicates that j a temperature variations that are observed over the 
the 1 1 / 12 is clearly the only period with such obvious con tinental regions. When an MSU undersamples this 
periodicity in the difference time series. area re i at ive to other areas, there is a potential for dra- 

Upon closer study we found that for a given latitude, tnatically increased noise in the difference between two 
the NOAA-12 MSU generates warmer than expected observing satellites. The spikiness of the time senes in 
temperatures in local winter and colder than expected 2 i s strongly related to this. Rather than averaging 

in local summer. In other words, the annual range of eac j^ zonal anomaly across a few days to reduce the 
temperatures in the NOAA-12 MSU was only about noise we shall use lhe tec h n ique of median hltenng in 
98% of the other MSUs. This observation might be w hi c h the outliers or spikes are essentially eliminated 
related to the fact that the dynamic range of the N0.4A- f rom the dataset 

12 MSU channel 2 was reduced due to much reduced For a given latitude and given day, the median is 
gain ( Mo 1994). determined from the current value ranked with n days 

To account for this effect, that is, adjust the dynamic previous and following. The median window* is then 
range of the NOAA-12 and possibly other MSUs, we 2 n + 1 days in length. Data will be set to missing if the 
decompose, latitude by latitude, the annual cycle of 2 n + 1 window finds less than n + 1 days of data. In 
the difference time series into Fourier components for pjg 5 we s h ow t he result of noise reduction to the daily 
those overlap cases longer than 1 year(NOAA 6/7, 6/ global-mean anomaly for no median filtering through 
9. 10/11, 11/12). Once these Fourier components of an j \_day window (n = 5) once the annual harmonics 
the residual are removed, the time senes is recon- h ave b een removed as described above. We notice a 
structed as before. The results for daily and 30-day substantial reduction in daily VjV for n~ 1 and 2, with 
means are shown in Fig. 4 for various Fourier-com- | esse r reductions thereafter. However, reducing .V is 
ponent wavenumbers removed (e.g., “5" indicates re- only one component of the goal, the other being a pres- 
moval of components 1 to 5). ervation of the signal magnitude S. For the same cat- 

The removal of Fourier components 0, 1, 1-3, 1-5, egories of median filtering, in Fig. 6 we display the 
1-7, and 1-9, based on the difference time series, pro- values of V 5 where reductions are also apparent in the 
duces S/X values for daily anomalies of 8.7, 9.2, 9.8, <j a j|y variance. A portion of the reduction in both A' 
9.9, 10.0, and 10.0, respectively. For 30-day averages, aiK j $ | s d ue to the removal of the errors in the times 
these values are 26.2, 30.6, 43.3. 45.2, 45.1, and 45. 1 . 

This adjustment impacts the 30-day S/ N substantially 
more than the daily, reducing )[X by over 20% in the 
case of waves 1-5. The daily VA r is also reduced, but 
less than 10% because its fluctuations occur more on 
a day-to-day basis. There is no significant reduction of 
X beyond wave l -5 removal. Given this magnitude of 9 
noise reduction we believe this procedure should be 
applied ( up to annual harmonic 5 ) to the data supplied 
to the users. The process of removing Fourier com- 
ponents of the annual period does not affect the overall 
trend of the time series at the 0.002°C per decade level 

and has its primary impact on NOAA-12. Fig. 6 - As above tor Vs. 
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Fig. 7. As in Fig. 2 for the SOAA-ll and SOAA-l 2 overlap com- 
parison after removal of five Founer components of the annual period 
and inclusion of 5-point median filtering. Note the trend in the dif- 
ference time senes. 


series. However, the natural smoothing characteristic 
of median filtering is also a factor. 

We believe it is desirable to include median filtering 
of the latitudinal anomalies because of the obvious re- 
duction of N. The selection of the width of the window 
is somewhat arbitrary since the level of signal smooth- 
ing one might accept varies according to research re- 
quirements. We Jiave chosen n = 2 for a width of 5, 
which reduces fS by only 5.6% while reducing VAfby 
33%. The 5/A r ratios of daily global anomalies for the 
median filter windows shown in Figs. 5 and 6 (width: 
2 n+ 1 = 0. 3, 5, 7, 9, and 1 1 ) are 9.9, 18.1, 26.2, 3 1 .8, 
38.1, and 41.1, respectively. Reductions in VS for the 
same sequence of values are 0%, 3.8%, 5.6%, 7.0%, 
8.4%, and 9.8%. The S/N ratios for 30-day averages 
were 45 ± 1.5 with no discemable relationship between 
window width and S/N. 

A median filter of width 5 would not be expected to 
replace the center or current value with the calculated 
median one time out of five (20%) for a random dis- 
tribution. In other words, the time-centered value 
would be the actual median 20% of the time. For pro- 
cesses. which have discemable variations on time spans 
longer than 5 days, the percentage of nonreplacement 
would be greater as the current value would have a 
higher probability of being the actual median value if 
there were a trend over the 5-day period. 

Analysis of the 5-day median replacement rate of 
the latitudinal anomalies indicates that about 35% are 
not replaced at high latitudes, with that percentage de- 
creasing to slightly more than 20% at the equator. This 
implies that tropical variations of zonal-mean temper- 
ature are so slow as to represent essentially random 
distributions on 5-day intervals. Concerning the almost 
80% that are replaced in The Tropics, the difference 
between the current and median (i.e., replacement) 
value produces a o of about 0.2°C. In higher latitudes, 
the 65% current and median values differ with a a of 
about 0.4°C. The effect of the latitudinal median fitter 
is rather slight on the global mean as compensation 
between latitude bands occurs so that the VS is reduced 
by only 5.6%. 

The reduction of Va from a combination of har- 


monic removal and median filtering is over 45% for 
daily observations and over 20% for the 30-day aver- 
ages. Figure 7 shows the same overlap of 11/12 as 
shown in Fig. 2 after these two procedures are applied. 
For this particular overlap example, the S/N for daily 
and 30-day averages are now 18.4 and 36.8, respec- 
tively, up from 5.8 and 13.6. Using the various window 
widths, we calculated the trend of the entire time senes 
and found variations of the trend are less than 0.002°C 
per decade from that of the selected {In A- l =5) filter. 

d. The issue of trends in the di fferences 

What is now apparent in the 11/12 overlap time 
series of differences ( Fig. 7 ) is the disturbing presence 
of an obvious trend. In actuality, the trend began in 
1990 when the temperatures of the 10/11 overlap be- 
gan to diverge. During the first 2 years of the 10/11 
overlap, there was no trend in the difference time senes 
greater than 0.01 °C per year (Christy and Goodridge 
1994). Since 1990, NOAA-11 observes a global tem- 
perature trend, which is 0.03°-0.04°C per year warmer 
than that of post- 1 989 NOAA-10 and all of A OAA-1 2. 

Because N0.4A-1 1 is in disagreement with two other 
satellites, our first assumption is that NOAA-1 1 is the 
source of the error. This suspicion is substantiated by 
NOAA-1 Ts orbital drift from the initial equatorial 
crossing time of 0130/ 1330 to 0430/ 1630 as ol Feb- 
ruary 1994 (Fig. 8). The question now becomes how 
does the global-tropospheric temperature relate to the 
changing time of observation ( i.e., diurnal cycle ) The 
troposphere is generally warmer at 1630 than at 1 330 
but is this warming entirely offset by the cooling that 
would be expected in the nighttime pass as it dnfted 
from 0130 to 0430? 

To provide further insight into this issue we calcu- 
lated the trend latitude-by-latitude for the difference 
between NOAA-11 and the two morning orbiters. 
NO A A - / 0 and - / 2 ( Fig. 9). We see that the subtropical 
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Fig. 8. Ascending pass crossing time at the equator tor SOAA-l ! 
from September 1988 to February 1994. 
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Latitude 


FIG. 9. Relative trend (°C per year) between A OAA-ll versus 
\OAA-IO and A OAA-I2 determined from 1 January 1990 to 28 
February 1994 as a function of latitude. 

latitudes show slight positive to negative trends for the 
relative difference, while at all other latitudes the sign 
is positive. The notion that the MSU on NOAA-I 1 is 
contributing a spurious warming trend related to its 
orbit time decay seems confirmed by this figure. In 
those regions with strong convection, in which surface 
heat is transferred to the free atmosphere throughout 
the afternoon, the warming measured from 1330 to 
1630 would be greater than cooling during the stable 
period from 0130 to 0430. The Northern Hemisphere, 
with its greater proportion of land surface relative to 
the Southern Hemisphere, seems to have a greater pro- 
pensity for vertical mixing and, therefore, generally 
supports more positive trends during the 01 30/ 1 330- 
0430/1630 drift. Even over those oceanic regions, 
which are subject to frequent precipitating systems ( i.e., 
tropics and midlatitudes), it has been shown that the 
tropospheric temperature peaks at about 1600 local 
time (Randall et al. 1991). 

However, in the subtropical latitudes there is con- 
siderable subsidence and, thus, decoupling between the 
boundary layer (most often a marine boundary layer 
in the Southern Hemisphere ) and the free atmosphere 
(Trenberth et al. 1992). The diurnal cycle in the tro- 
posphere would more closely follow the solar cycle, 
peaking in the early afternoon and cooling until dawn. 
Tropospheric temperatures at 0430 would certainly be 
cooler than those at 0130, while in the afternoon, the 
1630 observation would most likely be equal to or 
cooler than the 1330 observation. 

The latitude bands in the 32.5°-40°N region have 
a peculiar character in that they are somewhat affected 
by the typical subtropical atmospheric processes and 
also by pronounced surface emissivity effects. The 
highlands from Iran eastward to China penetrate the 
MSU weighting function, “shining through 1 ' the at- 
mosphere to some extent. The diurnal cycle of the MSU 
2R in these areas is very closely tied to the surface skin 
temperature, so that the MSU brightness temperature 
reaches its maximum earlier (near solar noon) than 
does the physical temperature of the air above. By 
viewing at 1630 rather than 1330, XOAA-ll observes 
a cooler brightness temperature due to the cooler skin 


temperature even though the atmosphere there may 
be warmer. The removal of a negative trend in those 
bands is, therefore, necessary. 

The effect of the satellite orbit decay would be greater 
during summer, when convection (dry and moist) is 
more vigorous, than during winter. As such, one would 
expect an annual cycle in the time-drift effect on the 
difference time series. This is not seen in our data at 
this point because such oscillations, which are anchored 
to the annual period, have already been removed ( sec- 
tion 4b). In fact, that these harmonics were present to 
begin with seems related to this orbit-drift problem as 
well as instrument differences. 

The trends of Fig. 9 are removed from the NOAA- 
II data, beginning with 1 January 1990. This process 
further reduces the noise of the last two overlap 
pairs and. thus, the overall noise levels. The final ViVand 
S/N for all overlapping comparisons are 0.0348 K 
(0.0167 K) and 29.5 (63.4) for daily ( 30 day averages ) 
respectively. The reductions of noise for the procedures 
described in this paper are summarized in Fig. 10. 

Our definition of noise, as the variance of the dif- 
ference time series divided by two, is not a traditional 
quantity but did suit our purposes for the visual com- 
parison of signal (the sum divided by two) and noise 
for intercompared satellite data. The standard error of 
measurement, as another representation of noise, re- 
quires that our Vfr be multiplied by 1.414 (i.e., 
2/V2). We, therefore, compute the daily (30 day av- 
erage) value of the standard error of measurement as 
0.0492 K (0.0236 K). 

In Fig. 1 1 we show the 11/12 overlap in which now 
there is no trend in the difference time series. We note 
here that the latitudinal trends removed in NOAA-I 1 
are linear functions of best-fit. We assumed that the 
actual warming/cooling rates due to time drift are lin- 
ear over the 3-hour drift observed in NOAA-I / As it 
continues to drift, we expect that the instrument will 
eventually observe temperatures at their local maxi- 
mum ( globally speaking) of the diurnal cycle and then 
the temperatures will begin to decline. Since NOAA- 
1 1 is degrading at 6 minutes per month, we suspect 


0.07 - 



Fig. 10. Solid (open) columns represent VN for daily (30-day av- 
erages) for A: unaltered data for 20 k observations. B. as in A with 
removal of five Fourier components. C: as in B with 5-point median 
tiltenng. and D: as in C after removal of NOAA-I l trend. 
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NOAA-1 1 and 12 
Daily Global Average Anomaly 
September 1991 to February 1994 
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Fig. 1 1. As in Fig. 7 for the AOcL-J-7 / and - 12 overlap after latitudinal 
trends ot' NOAA-1 1 are removed. 


that the usefulness of the linear adjustments made here 
may have already reached their limit. 

The orbital drift in the afternoon satellites is not 
detectible in the temperature data until the third year 
of orbit, as in the case of NOAA-1 1. For TIROS-N, 
there was a sudden drop in the MSU data rate after 
only 14 months (January 1980), so drifting had no 
impact. The oscillator on NOAA-9 ' s MSU failed after 
only 27 months, rendenng the instrument inoperative 
for channel 2. We should note that other instruments 
remained in operation on NOAA-9 and researchers us- 
ing those data will certainly have to deal with the prob- 
lem of drift. The MSU on NOAA-7 was in operation 
for 3.5 years, and, as such, there may be a slight warm- 
ing effect late in its term of service (about 1984, Fig. 
1 ) when it was the only MSU in operation. We should 
note, however, that the subsequent MSU on NOAA-9 
was calibrated directly against NOAA-6, and so did not 
require NOAA-7 data. The long time series of NOAA- 
6 was rigorously checked by radiosonde data from mid- 
1979 through 1986 and found to have no temperature 
drift over the 7.5 years ( Spencer etal. 1990). Therefore, 


the MSU time series after NOAA-7 is independent from 
any temperature drift that might have occurred be- 
tween the first shutdown of NOAA-6 ( Apnl 1983 ) and 
its restart in April 1985. 

5. Conclusions 

Careful examination of the daily global-mean 
anomalies of the lower-troposphenc temperature 
through the comparison of simultaneously operating 
satellites has allowed a direct determination of signal 
and noise of the data. From this, we have determined 
which characteristics are important in noise reduction. 
In the case of missing data, there was little impact found 
between the 19 000 and 22 000 per day data rates, so 
the value of 20 000 was retained. 

There were significant impacts due to residual annual 
cycle differences between the various instruments and 
the lack of latitudinal /temporal filtering. By removing 
the residual annual cycle to 5 Fourier components and 
applying a median filter of width 5 days tor each lati- 
tude, the daily and monthly noise were reduced by 
over 40% and 25%, respectively, while leaving the signal 
with over 90% of its original variance. 

The time drift of NOAA-1 1 was shown to produce, 
since 1990, a spurious warming trend in the time series 
as it observed the atmosphere at warmer local times 
(globally averaged) during its degradation. Once cal- 
culated and removed, the resultant total time senes 
produces a decadal trend for the 15+ year time scries 
of-0.058°C (Fig. 12). This is 0.03 °C per decade o^ier 
than calculated under the assumption that \(> \ i / 1 
was in a stable orbit. Also in Fig. 12 is the time scries 
of the differences between previously reported global 
temperatures and those using the adjustments applied 
from this current research. Monthly differences are cr\ 
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Fig. 12. Time senes of monthly global anomalies (thick line) relative to 1982-91 annual cycle alter 
removal of effects noted in this paper. Difference (thin line, offset by -0.6°C for clanty) relative lo MSU 
2R values reported prior to adjustments given in this paper. The time senes here includes March 1994. 
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minor (within 0.03°C) until the last 3 years when the 
impact of the NOAA-1I drift is apparent. 

In a recent report, Christy and McNider ( 1994) es- 
timated the magnitude of the global MSU 2R temper- 
ature trend once the effects of volcanoes and El Ninos 
were removed. The goal was to isolate a possible 
warming signal due to the enhanced greenhouse effect. 
Their resultant 15-year trend was about -f0.09°C per 
decade. Using the MSU 2R data, as adjusted in this 
study, the volcanic impact was reevaluated because the 
impact of Pinatubo (June 1991) occurred entirely 
within the period, which was affected by the NOAA- 
11 orbital drift. The global temperatures since 1991, 
now cooler than previously thought, imply that Pina- 
tubo had an even stronger impact on the global tem- 
perature than estimated in Christy and McNider. 
However, after reevaluating the volcanic effect on the 
now-cooler post- 1990 global temperatures and then 
removing that effect we find the "global warming" 
trend remains at about +0.09°C per decade as of July 
1994. 

An opportunity is now available to assess the diurnal 
cycle of temperature of the global- (and zonal) free 
atmosphere due to the inadvertent drift of NOAA-1 1. 
With the morning satellites well anchored to the ter- 
minator orbit, we can compare global and zonal tem- 
peratures in NOAA-1 1 with the morning orbiter data 
throughout the 3-hour drift as was done in this paper. 
From this, we may determine average changes between 
the drifting observations and those of the morning or- 
biters to assess part of the diurnal cycle. Investigations 
concerning the diurnal cycle, as well as correcting MSU 
channel 4 (lower stratosphere), will be addressed in 
future work. 
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